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Abstract 

QP ■ Different from previous models based on scatter theory and random matrix theory, a new in- 
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terpretation of the observed log-normal type time-headway distribution of vehicles is presented in 
this paper. Inspired by the well known Galton Board, this model views driver's velocity adjusting 
process similar to the dynamics of a particle falling down a board and being deviated at decision 
points. A new car-following model based on this idea is proposed to reproduce the observed traffic 



?H i flow phenomena. The agreement between the empirical observations and the simulation results 

j^ ' suggests the soundness of this new approach. 
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I. INTRODUCTION 

To explain and reproduce the complex phenomena of road traffic, the dynamics of traffic 
lows are often described on A^ strongly-linked particles (vehicles) under fluctuations [l|], 2], 
31 . Since the eroverninEf interaction forces or potentials cannot be directly measured, the 

n n 

statistical distributions of particles are often investigated instead [J], [5|. 

Among different statistical features, the distributions of space-gaps/time-headways be- 
tween these particles (vehicles) received consistent interests from various viewpoints. In 
transportation engineering, the corresponding investigation data are helpful not only in un- 
derstanding the microscopic-level driving behaviors but also in estimating the macroscopic- 



level roadway capacities [6j 



0,B, 



9|. Numerous distribution models have been developed 



over the past 50 years to directly fit the empirical data. In some recent studies ll|, 12|, 13| . 
|l4| . some theoretical models were presented from various physical perspectives (e.g. scatter 
theory and random matrix theory) to explain why we can observe similar distributions even 
or different phases (i.e. Kerner's the free-flow, synchronized flow, and moving jam phases 



Differently in this paper, another interpretation inspired by the well-known Galton board 
is presented. A new car-following model based on this interpretation is also proposed to 
reproduce the observed traffic flow phenomena. The research purposes are twofold here: 
1) Previous physical interpretations focus on the steady-state macroscopic-level statistics; 
while this method provides a microscopic-level dynamic explanation, which can also be 
used to simulation the transient-state statistics of inter-arrival and inter-departure vehicle 
queuing interactions. 2) It is interesting that as a slightly modified extension of the clas- 
sic car-following model, this new model can easily reproduce the observed time-headway 
distributions, which had been neglected by many previous approaches. 

II. THE LOG-NORMAL DISTRIBUTION MODEL OF TIME-HEADWAYS 

Based on several recent studies [5|], |9|, [lO|, [ll[, [ij], we believe that log-normal dis- 
tribution model is a simple yet effective model comparing to other ones, i.e. distribution 
models according to random- matrix theory |5|, [ij]. For example. Fig. [T 
ison results for the empirical, super-statistical on random matrix theory 



shows the compar- 

n 

14l | and log-normal 



distributions P{t), where r = th/{th), th denotes the sampled time headways. It is clear 
that log-normal distribution model leads to smaller fitting errors. 
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FIG. 1: Probability density P{t) for scaled empirical time headways r between successive cars 
in traffic flow. For comparison, the empirical data and super-statistical distributions on random 
matrix theory are excerpted from [5|]. 



Moreover, log-normal distribution model can also be used to depict the distributions of 
time headways in transient-state flow. For instance, the departure headways (which are 
usually defined as the times that elapse between consecutive vehicles when vehicles in a 
queue start crossing the stop line) are shown to follow log- normal distributions in [l9|. 

Noticing that the random matrix theory gives a soundable explanation for its corre- 
sponding distribution model, an interesting question naturally arises as: "Can we provide a 
physical interpretation of such log-normal distributions?" After comparing with the other 
dynamic processes that yield log- normal distribution [l5|, [16|, 17|, 18|, we think the out- 
comes of such distributions can be explained as follows. 

In car-following process, the driver of the following vehicle will adjust his/her velocity 
from time to time to track the leading vehicle and meanwhile keep a safe distance between the 
leading vehicle and him/her. Because the leading vehicle's movement is often unpredictable 



(at least not fully predictable), the accelerating and braking action of the driver is often 
overdue. 

This persistent velocity adjusting process is somewhat like the process of the particles 
falling down a board and being deviated at decision points (the tips of the triangular ob- 
stacles) either left or right with equal probability, see Fig. [2J If the deviation of the particle 
from one row to the next is a random additive process with possible values +c and — c, 
the normal distribution will be created by the board, which reflects the cumulative additive 
effects of the sequence of decision points. But if the deviation of the particle from one row to 
the next is a random multiplicative process with possible values -c! and /c', the log-normal 
distribution will be generated. 

It is obvious that drivers tend to take very careful acceleration when the spacing between 
the two vehicles is small and do not want to speed up at once after braking. Thus, the de- 
viation of the particle to the left is equivalent to the following vehicle's accelerating action, 
since the relative speed decrease slowly at this period; while the deviation to the right is 
equivalent to the following vehicle's braking action, since the relative speed increase quickly 
at this period. From this view point, the car- following process is similar to the particles 
falling in the log-normal type Galton board. To verify this conjecture, a microscopic sim- 
ulation model inspired by log-normal type Galton board is proposed in the next section to 
reproduce the observed phenomena. 

III. THE NEW CAR-FOLLOWING MODEL INSPIRED BY GALTON BOARD 

In this section, it will be interesting to find that a car-following model incorporating 
Galton Board mechanism can be easily derived from a famous classic car-following model 
proposed in later 1950s [20|, |2l|]. The proposed model mainly depicts three states of a 
vehicle: 1) stopped, 2) starting-up, and 3) driving/braking states. More precisely, they can 
be described as follows. 

Suppose Vi{t) and Xi{t) are the velocity and position of the iih vehicle (follower) at time t, 
respectively. Similarly, Vi-i{t) and Xj_i(t) are the velocity and position of the i — 1th vehicle 
(leader) at time t. gi{t) = Xi{t) — Xi-i{t) — U denotes the gap between the two vehicles at 
time t. T > is the simulation time span. 





(b) 

FIG. 2: Diagram of two Galton boards yielding normal (a) and log-normal (b) distributions re- 
spectively. If the tip of a triangle is at distance x from the left edge of the board, triangle tips to 
the right and to the left below it are placed at a; -|- c and x — c for the normal distribution, and 
X ■ c' and x/c' for the log-normal, where c and c' are constants. 



1) If Vi{t) = 0, that is, the ith vehicle is fully stopped at time t; it will stay stopped in 
the next time interval T and then enter the starting state since time {t + T). 

2) Else if the ith vehicle is in the starting-up state at time t, it will check its coming 
actions from the following three choices: 

2. a) If Qiit) < Gmini whcrc Gmin IS the minimum stop distance. The zth vehicle 
will go back to stop state by letting Vi{t + T) = without doing anything else (that's, 
Xi{t + T) = Xiit). 

2.b) Else if Vi{t) < v starting, wheie v starting IS the limit of the staring-up velocity; the 



vehicle is still in the starting-up state at time (t + T). Thus, it will accelerate with the 
staring- up rate agtarUng in the next time interval T till time (t + T). 

2.c) Otherwise, the vehicle will switch to the driving/braking state since time t and 
re-determine its action according to the updating rules set for the driving/braking state at 
time t. 

3) Otherwise, the vehicle is running in the driving/braking state. Actually, the driv- 
ing/braking state contains three modes: 3. a) free-driving mode, 3.b) braking model and 3.c) 
following mode (see Fig. [3]), which can be represented in details as: 



Free-Driving Mode 




Vi(t) - Vi_i(t) 



FIG. 3: Diagram of the three modes (regions) in the gap vs. relative-speed space. Comparing to 
Fritzsche 1994 model used in PARAMICS and Wiedemann Psychophysical model used in VISSIM, 
our model has fewer modes. 

3. a) If gi(t) > Gmax, where Gmax is the maximum coupling distance; the ith vehicle is 
in free-driving mode. It will try to approach the highest velocity as 



Vi{t + T) = mm{vi{t) + a+„^. x T, Vmax} 



(1) 



where amax is the maximum allowable acceleration rate. For simplicity, we choose T = 1 
second here. 

3.b) Else if fj(t) — fj_i(t) > [gi{t) —G]/H, there exists a risk to collide, and ith vehicle 
is in braking mode. It will try to stop as quickly as possible 



Vi{t + T) = v,{t) - D xT 



(2) 



where G, H and D are positive constants denoting the minimum safety gap, deceleration 
time and braking decelerating rate, respectively. 

3.c) Otherwise, the ith vehicle is in the following mode. We have the following two-step 
updating rules as 



Vi{t + T) 



m.ax{(3vi{t) ,^_rp^_^\ i^j^^i. , 0}, with probability p to decelerate 



' Xi(t-T)-x,_i{t-T)-li 

Xi{t)-Xi-l{t)-li 

'xi{t-T)-x,_i{t-T)~~k^ 



mm{^Vi{t) ^ f^^lp^_l' ]f^_r)^j , Vmax}, with probability (1 - p) to accelerate 

(3) 



. , ^, . max{vi{t),v,_i{t)-a^^^^^xT}, if decelerate 
Vi{t + T) = ^ (4) 

rmn{vi{t) , Vi-i{t) + a^ax ^ ^l^ if accelerate 
where U denotes the length of the ith. vehicle. Vmax is the maximum allowable velocity, a^^^, 
and a~^3, are the maximum allowable car-following ac/deceleration rates. The max, min 
functions are added to guarantee that the velocity and ac/decelerating rates are within the 
limits. /5 is a positive constant to be determined according to investigation data. 

Eq.([3])-(j4]) characterize the simulated time headways to follow the log-normal distribu- 
tions. A brief explanation of this mechanism is provided in the Appendix. It is assumed 
that < /5 < 1 together with < j> < 0.5, since drivers tend to keep a close spacing. We can 
see that Eq.(l3])-(j4l) naturally incorporate the frequently-used randomly- slow- down features 
of driving behaviors. 

4) Finally, the position of vehicle is updated as 

Xi{t + T)=Xi{t)+Vi{t + T)xT (5) 

To illustrate the effectiveness of the proposed model, a specialized form is used to repro- 
duce the complex phenomena of freeway traffic flows. Particularly, we have all the vehicles 
have the same length as U = 4m and T = Is. Other parameters are set as v starting = 8m/s, 
astarting = 4m/s^ Gmin = 0.5m, Gmax = 52.5m, Vmax = 30m/s, 0+^3, = 4m/s, a;;„3, = 8m/s, 
P = 0.93, p = 0.3. G = 0.5m, H = 7s, D = Sm/s^. 

Fig. m shows the corresponding distributions of the simulated time headways within dif- 
ferent velocity ranges. The log-normal distribution model passes the Kolmogorov-Smirnov 



(K-S)hypothesis test. This proves that the new model yields the log-normal type time- 
headway distribution as observed; while many previous models yield symmetric distribu- 
tions, i.e. Fig. 7 of |23|. Moreover, as the velocity of the vehicle increases, the mean time 
headways of such log-normal type distributions will approach the saturation headway. This 
fits the observations [5(], [2J], too. 
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FIG. 4: Probability density for (a) scaled simulated time headways r and (b) un-scaled simulated 
time headways th between successive cars in traffic flow. 



Fig. [5] shows the corresponding distributions of the simulated space gaps under different 



traffic pressures. As the traffic flow density (occupancy) increases, the distribution of space 
gaps will become more concentrated, which is in accordance with the practical observations 
(see the un-scaled distributions shown in Fig. 3a in |13[ for comparison). Moreover, the 
shape of the un-scaled simulated gap distributions look similar to the empirical data and 
also the scaled distributions predicted by using random matrix theory (see that shown in 
Fig. 1 in [ij] for comparison). 
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FIG. 5: Probability density for (a) scaled simulated space gaps l = gs/ids) and (b) un-scaled 
simulated space gaps gs between successive cars in traffic flow. 



Fig. [6] show the {k, Vg) and {k, q) diagrams for the proposed model, obtained by local and 
global measurements. These macroscopic traffic stream characteristics are measured accord- 
ing to 



22| . Here, Kg denotes the length of the closed single-lane system, which is 27000m. 
It can be found that the local measurements in Fig. [H] also discriminate Kerner's three-phase 
flows: the free-flow regime contains only a few data points on a line; the synchronized regime 
is formulated by widely scatters of the data points; and jammed regimes contains the data 
points corresponding to Kerner's line J. 

Fig. [7] shows three typical complex phenomena reproduced by using the proposed model: 
(a) the oscillation in traffic flow; (b) the wide moving jam; and (c) the stop-and-go waves. 
This proves that the proposed model is capable to describe the other complex dynamic 
features of road traffic flow, too. 

APPENDIX: EXPLANATION OF THE LOG-NORMALITY 

According to the Galton board model, if the variation dynamics of time headway should 
be depicted by the following equation without losing any generality. 

, ,,. J Pth{t-T), with probability (1 - p) 
^th{t — T), with probability p 
Notice that be definition, we roughly have 



tnit) - -f%. (A.2) 

v{t + T) 



we can then have 



9s{t) J /^^4?r' with probability (1-p) 



'"^^^^^ 1 y-"^^^ with probability p 



(A.3) 



Considering the pre-determined limits of the velocity and ac/deceleration rates, we can 
directly derive Eq.(jll) from Eq.( lA.3l) then. 

ACKNOWLEDGMENTS 

This work was supported in part by National Basic Research Program of China (973 
Project) 2006CB705506, Hi- Tech Research and Development Program of China (863 

10 



35 



Globle Measurement 

o Local Measurement 




500 1000 1500 2000 2500 3000 3500 4000 

Vehicle Number 



(a) 




500 



1000 1500 2000 2500 3000 3500 4000 
Vehicle Number 



(b) 

FIG. 6: (a) The {k, Vg) diagram for the proposed model, obtained by local and global measurements, 
(b) The {k,q) diagram for the proposed model, obtained by local and global measurements. 
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